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Abstract 

This short article describes flow parameters, numerical method, and 
animations of the fluid dynamics video "LES of an Inclined Jet into 
a Supersonic Turbulent Crossflow" ( high-resolution] and |low-resolution| 
video). We performed large-eddy simulation with the sub-grid scale 
(LES-SGS) stretched-vortex model of momentum and scalar transport 
to study the gas-dynamics interactions of a helium inclined round jet 
into a supersonic (M = 3.6) turbulent (Reg = 13 x 10 3 ) air flow over 
a flat surface. The video shows the temporal development of Mach- 
number and magnitude of density-gradient in the mid-span plane, and 
isosurface of helium mass-fraction and A2 (vortical structures). The 
identified vortical structures are sheets, tilted tubes, and discontinuous 
rings. The vortical structures are shown to be well correlated in space 
and time with helium mass-fraction isosurface (Yfje = 0.25). 

1 Flow parameters 

Figure [T] shows the flow schematic. Helium is injected through an inclined 
round jet into a supersonic turbulent air flow over a flat surface. In the 
present investigation, the jet axis forms a 30° angle with the streamwise 
direction of the air flow. The flow parameters of air and helium are reported 
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in Table [l] The jet diameter, d, is 3.23x10 3 m, and the boundary layer 
thickness, 5, of the air flow is 2xl0~ 2 m as in the experimental study of 



Maddalena, Campioli & Schetz (2006). The air free-stream Mach number is 



3.6, the jet Mach number is 1.0, and the jet to free-stream momentum ratio, 
q = {pU 2 )-J (pU 2 ) e , is 1.75. The Reynolds number of the air flow based on 
the momentum thickness is Reg = \J G Qjv v = 13 x 10 3 (Res = C/ e 5/f w = 
113 x 10 3 ), where U e is the free-stream air velocity and z/ w is the kinematic 
viscosity of air computed at the wall for adiabatic wall conditions. 
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Figure 1 : Flow schematic. 
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Table 1 : Flow parameters. 

2 Numerical method 

Large-eddy simulation was performed using the sub-grid scale (LES-SGS) 
stretched-vortex model of momentum and scalar transport developed by 



Pullin and co-workers ( 


Misra & Pullin 


1997 


Voelkl et al. 


2000||Pullin 


2000 


Kosovic et al. 


2002 


) . A hybrid numerical approach ( 


Hill&Pullin[|2004 


Pan- 



tano et al. 2007) with low numerical dissipation that uses tuned centered 



finite differences (TCD) in smooth flow regions, and a weighted essentially 



non-oscillatory (WENO) scheme (Liu et al. , 1994 Jiang & Shu, 1996 ) around 
discontinuities and ghost-fluid boundaries is employed. The level-set ap- 



proach with the ghost-fluid method (Fedkiw et al. , 1999 ) is used to treat the 



complex boundary (Fig. [TJ where no-slip and adiabatic boundary conditions 
are applied. The framework for block-structured Adaptive Mesh Refine- 



ment in Object-oriented C++ (AMROC) by Deiterding (2003) was adopted. 
The computational domain is a parallelepiped with sides —0.025 m < x < 
0.0774 m, -0.006 m < y < 0.0452 m, and -0.0256 m < z < 0.0256 m (i.e., 
-7.74 < x/d < 24, -1.86 < y/d < 14, and -7.92 < z/d < 7.92), in the 
streamwise, wall-normal and spanwise directions, respectively. The flat-wall 
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boundary is the y = plane and the center of the jet-exit plane is at the 
axes origin (Fig. [TJ. The basic mesh is 256 x 128 x 128 = 4.2 x 10 6 grid 
points with a base mesh spacing of Axb = 4 x 10 -4 m. The dynamic mesh 
refinement (AMR=2) adds one level of finer mesh on the coarse mesh, in- 
creasing the effective total number of mesh points to approximately 16 x 10 6 
with a fine-mesh spacing of Axf = 2 x 10~ 4 m. The use of AMR=2 saves 
about half grid points with respect to the fully refined mesh of 33 x 10 6 grid 
points. The Kolmogorov length-scale, rj, of the shear layer created between 
the high-speed helium jet and low-speed air stream is of about 1.5 fim, i.e., rj 
is of 0(100) times smaller than the fine-mesh spacing. The SGS-TKE is less 
than 20% of the total TKE in most of the flow field . Thus, the resolution 



criterion for a sufficiently resolved LES-SGS ( |Pope 2004 1 is satisfied. 
2.1 Inflow conditions 

The transition and spatial development of the helium jet were found to 



be strongly dependent on the inflow conditions of the crossflow (Ferrante 



et al. 2009-1511). These results indicate that correct turbulent inflow con- 
ditions are necessary to predict the main flow characteristics, dispersion 
and mixing of a gaseous jet in a supersonic turbulent crossflow. A method- 
ology for the generation of synthetic turbulent inflow conditions for LES 
of spatially developing, supersonic, turbulent wall-bounded flows has been 



developed by Ferrante et al. (2010). A brief description is here given. In- 



flow conditions are computed as the sum of zero-pressure gradient, com- 
pressible turbulent-boundary-layer (ZPG-CTBL) mean profiles and turbu- 
lence fluctuations. First, the friction velocity, u T , and velocity profile (von 



Karman 1930) of the incompressible ZPG-TBL are computed by prescrib- 



ing the Reynolds number based on boundary-layer thickness, Res (Ferrante 



& Elghobashi 2004). The resulting velocity profile is then transformed into 
the velocity profile of the ZPG-CTBL at M e = 3.6 according to van Driest 
(Ivan Driestl 119511 ISmits & Dussaugel 120061 iKnightl 12006-4981). Last, the 



temperature profile of the ZPG-CTBL is computed using the Walz formula 



(Walz 1969) assuming an adiabatic wall. Inflow turbulence fluctuations are 



generated by modifying the methodology of (Ferrante & Elghobashi 12004) 
to supersonic flows at high-Reynolds number. A model spectrum (Pope 



2000 )of turbulence kinetic energy, E(k), and the Reynolds stresses, (uiUj)^, 
are prescribed at the inflow plane. The Reynolds stresses at M e = 3.6 are 
obtained by scaling the incompressible Reynolds stresses of a ZPG-TBL at 
Re e =2900 obtained via DNS ( |Ferrante & Elghobashi} |2005[ ) with the the- 
oretical local profile p/p w . Such scaling is justified by the DNS results of 



Guarini et al. (2000) 
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3 Animations 



The fluid dynamics video "LES of an Inclined Jet into a Supersonic Turbu- 
lent Crossflow" (high-resolution and l ow- resolution! video) shows five anima- 
tions (Al to A5)Q 

1. Mach-number, M, contours in the mid-span plane; 

2. contours of density-gradient magnitude, |Vp|, in the mid-span plane; 

3. isosurface of helium mass-fraction, Yn e = 0.25; 

4. isosurface of A2 (vortical structures Q 

5. overlapped isosurface of helium mass-fraction Ijie = 0.25 (yellow), and 
vortical structures (blue). 

The Mach-number contours (Al) show that boundary-layer turbulence/ 
bow-shock interaction results in shock-wave unsteadiness that affects the 
roll-up of the shear-layer (formed between air stream and helium jet) and, 
consequently, modulates in time the size and shape of the barrel-shock (black 
region near the jet exit). The contours of the density-gradient magnitude 
(A2) show the bow-shock and the shear-layer formed in between the air 
stream after the bow-shock and the expanded helium jet. Both contours 
show large-scale structures advected downstream. Jet unsteadiness, lateral 
and wall-normal helium dispersion, and three-dimensional structure of the 
helium jet are shown in the animation of Yfi e = 0.25 isosurface (A3). The 
vortical structures (isosurface of A2 shown in A4) are sheets near the jet exit 
where the shear formed between the air-stream and the helium jet is large. 
Downstream the jet exit, the vortical structures are mostly tilted tubes which 
sometimes look like discontinuous rings. In A5, the isosurface of Yn e = 0.25 
(yellow) mostly envelopes the isosurface of A2 (blue). The two isosurfaces, 
showing helium-jet puffs and vortical structures, are well correlated in space 
and time. The vortical structures look like muscles that move the isosur- 
face of helium mass-fraction, contributing to helium dispersion and helium 
mass-fraction convoluted isosurface. 



1 Each animation lasts about 23 s and shows 30 frames/s for a total of 696 frames. The 
temporal development of the flow is shown for about 10 -3 s. Thus, the flow is shown 
23,000 times slower than at its actual speed. 



2 The vortical structures are educed using the A2-method (Jeong & Hussain 19951, 
where A2 is defined as the second largest eigenvalue of the tensor (SikSkj + flikQkj), 
where Stj = {djUi + diUj)/2 is the strain rate tensor, and = {djUi — diUj)/2 is the 
rotation rate tensor. 
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